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Abstract. We develop a Hamilton- Jacobi theory for singular la- 
grangian systems using the Gotay-Nester-Hinds constraint algo- 
rithm. The procedure works even if the system has secondary 
constraints. 
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1. Introduction 



One of the most classical problems of theoretical mechanics is the 
study of constrained systems. Essentially, there are two different mean- 
ings to understand constrained systems. One refers to systems where 
we externally impose constraints allowing some particular motions (ex- 
ternal constraints). The second case is when the degeneracy of a la- 
grangian function imposes constraints on the phase space of the system 
(internal constraints). In this paper, we will restrict ourselves to this 
last situation. 

At a first step, when the lagrangian is singular, there appear con- 
straints restricting the admissible positions and velocities. Later on, 
the evolution of these initial constraints may produce new constraints. 

The theory of degenerate (or singular) lagrangian systems is relevant 
in Field theory, and just the quantization of these systems led to Dirac 
[8] to develop a wonderful theory of constraints, later geometrized by 
Gotay, Nester and Hinds [HI EH E51 E3]- Recently, M. Leok and col- 
laborators [T7] have studied degenerate lagrangians arising from truly 
mechanical systems, even in presence of additional nonholonomic con- 
straints (see also the paper by de Leon and Martin de Diego [21J). 

Another important topic in theoretical mechanics is the Hamilton- 
Jacobi theory which allows us to find solutions of a hamiltonian systems 
by means of solutions of a partial differential equation, the Hamilton- 
Jacobi equation. Conversely, we can treat to solve a PDE using the 
characteristic curves of a Hamiltonian system (see these two standard 
books [H |2] for a general view of the theory and some modern ap- 
proaches in terms of lagrangian submanifolds; see also [28] for a more 
classical view) . In [T5J [20] , we have successfully extended the classical 
Hamilton- Jacobi theory for nonholonomic systems, and in [191 122] for 
classical field theories. Therefore, it seems quite relevant to extend the 
Hamilton- Jacobi theory also for degenerate lagrangian systems, and 
this is just the goal of the present paper. 

Briefly, the standard formulation of the Hamilton- Jacobi problem is 
to find a function S(t, q A ) (called the principal function) such that 



If we put S(t, q A ) = W(q A ) — tE, where E is a constant, then W 
satisfies 



W is called the characteristic function. Equations (11. ip and ( II. 2p are 
indistinctly referred as the Hamilton- Jacobi equations. 
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There have been several attempts to develop a Hamilton- Jacobi the- 
ory for degenerate lagrangian system ([2U [251 E7J)- These procedures 
were based on the homogeneization of the given lagrangian, which leads 
to a new lagrangian system with null energy; then, it is possible to dis- 
cuss the Hamilton- Jacobi equation for the constraints themselves. The 
main problem is that, due to the integrability condition for the re- 
sultant partial differential equation, one can only consider first class 
constraints. Therefore, the treatment of the cases when second class 
constraints appear should be developed by ad hoc arguments (as in 
[27] . for instance). Thus, in [2J] and [25] the authors only discuss the 
case of primary constraints. 

Therefore, the Hamilton- Jacobi problem for degenerate lagrangian 
is far to be solved. 

Our procedure to develop a geometric Hamilton- Jacobi theory is 
strongly inspired in two main issues. The first one in the recent ap- 
proach to the Hamilton- Jacobi theory developed by Carinena et al [5j 
(see also [61 [26] , and [20J for the applications to nonholonomic me- 
chanics and field theory); and the second one, is the geometric theory 
of constraints due to Gotay and Nester [15j . 

Let us recall that given an almost regular lagrangian L : TQ — > M 
one can define a presymplectic system on Mi = FL(TQ) C T*Q, 
the primary constraint submanifold where U\ is the restriction of the 
canonical symplectic form on T*Q to M\, and FL : TQ — > T*Q is the 
Legendre transformation defined by L. The dynamics is obtained from 
the equation 

ix<^i = dhi, 

where hi G C°°(Mi) is the projection of the energy El G C°°(TQ). 

The above equation produces a sequence of submanifolds 

• • • My. <—}•••<-)■ Mi M\ <^-> T*Q 

and, eventualy, a final constraint submanifold Mf if the algorithm sta- 
bilize at some step. 

The strategy is to consider the projection of the constraint subman- 
ifolds provided by the constraint algorithm, so that we obtain new 
surjective submersions onto submanifolds of the given configuration 
manifold. This fact permits to connect a given solution of the final con- 
straint submanifolds Mf, with its projection onto Qf {jif : Mf — )■ Qf 
is the surjective submersion) using a section of 7Tj. 

The SODE problem is also discussed such that one can obtain the 
corresponding lagrangian picture. 

We also discuss the relation of the geometric Hamilton- Jacobi prob- 
lem with the Hamilton- Jacobi problem (in a traditional sense) for ar- 
bitrary extensions of h\ , in terms of first and second class primary and 
secondary constraints. Therefore, this work can be considered as the 
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natural extension to the Hamilton- Jacobi problem of the geometriza- 
tion by Gotay and Nester of the Dirac constraint algorithm. 

Several examples are discussed along the paper in order to illustrate 
the theory. 



2. Classical Hamilton- Jacobi theory (geometric version) 

The standard formulation of the Hamilton- Jacobi problem is to find 
a function S(t, q A ) (called the principal function) such that 

f +*/,§) = <», 

where h = h(q A ,pA) is the hamiltonian function of the system. If we 
put S{t, q A ) = W{q A ) - tE, where E is a constant, then W satisfies 

h{qA i w ] = E] (2 - 2) 

W is called the characteristic function. 

Equations (12. ip and (12. 2p are indistinctly referred as the Hamilton- 
Jacobi equation. 

Let Q be the configuration manifold, and T*Q its cotangent bundle 
equipped with the canonical symplectic form 

loq = dq A A dp A 

where (q A ) are coordinates in Q and (q A ,PA) are the induced ones in 
T*Q. In what follows, ttq : T*Q — > Q will denote the canonical 
projection. 

Let h : T*Q — y R a hamiltonian function and Xh the corresponding 
hamiltonian vector field, say 

ix h = dh. 

Therefore, the integral curves (q A (t),pA(t)) of X h satisfy the Hamilton 
equations: 

dq A dh dpA dh 



dt dp a ' dt dq A ' 
We can define also the Poisson bracket ot two functions. Given / and 
g real functions on T*Q, we define a new function {/, g} by 

{f,g} = u Q {X f ,X g ), 

where Xf and X g are the corresponding hamiltonian vector fields. 

The Poisson bracket gives us the evolution of observables, since given 
the hamiltonian h we have 

/ = X h (f) = ix h (iXfU Q ) = uj Q (X f ,X h ) = {/, h}, 
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and then we can rewrite the Hamilton equations as 

q A = {q A , h} 
Pa = {pA,h}. 

Let A be a closed 1-form on Q, say dX = 0; (then, locally A = dW). 

The following theorem gives us the relation of the Hamilton- Jacobi 
equation and the solutions of the Hamilton equations (see [HE]). 

Theorem 2.1. The following conditions are equivalent: 

(i) If o : I — >■ Q satisfies the equation 

dq A dh 
dt dp a ' 

then A o a is a solution of the Hamilton equations; 

(ii) d(h o A) = 



We can reinterpret Theorem 12.11 as follows (see [3 [T8j 120]). 
Define a vector field on Q: 

X% = Ttcq o X h o A 
The following diagram illustrates the construction of the vector field 

T*Q — ^ T(T*Q) 




Q 



h 



TQ 



Then the following conditions are equivalent: 
(i) If a : I — > Q satisfies the equation 

dq A dh 
dt dp a ' 

then A o o is a solution of the Hamilton equations; 
(i)' If a : I — » Q is an integral curve of X£, then Aocr is an integral 

curve of Xh] 
(i)" Xh and X£ are A-related, i.e. 

TX(X^) = X h o X 



Next, we have the following intrinsic version of Theorem 12.11 

Theorem 2.2. Let X be a closed 1-form on Q. Then the following 
conditions are equivalent: 

(i) X^ and Xh are X-related; 

(ii) d{h o A) = 
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Proof: In local coordinates, we have that 

dh d dh d 



X h 



and 
Then, 



dp a dq A dq A dp a 
A = X A (q) dq A 



dp a dq A ' 

dh d dh d\s d 



dp a dq A dp a dq A dps 
dh dh d\s s , A 



Since d\ = if and only if 

d\ A _ d\ B 
dq B dq A ' 

we have the equivalences between (i) and (ii). □ 
If 

A = X A (q) dq A 
then the Hamilton- Jacobi equation becomes 

h(q A ,X A (q B )) = const. 
If A = dW then we recover the classical formulation 

, / a 9W^ 

h{q , -t^x) = const. 

since 

dW 

3. The Hamilton- Jacobi theory in the lagrangian setting 

Let L : TQ — > R be a lagrangian function, that is, 

L = L(q A ,q A ) 

where (q A , q A ) denotes the induced coordinates on the tangent bundle 
TQ of the configuration manifold Q. In what follows, tq : TQ — > Q 
will denote the canonical projection. 

Let us denote by 

S = dq A ~ 9 



dq A 
and 
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the vertical endomorphism and the Liouville vector field on TQ (see 
[23J for intrinsic definitions). 

The Poincare-Cartan 2-form is defined by 

ujl = —dctL , «l = S*(dL) 
and the energy function 

E L = A(L) - L 
which in local coordinates read as 

«l = Pa dq A 
Ul = dq A dpA 
E L = q A p A ~ L(q, q) 

dL 

where ft a = tt-t stand for the generalized momenta. Here S* denotes 
oq A 

the adjoint operator of S. 

The lagrangian L is said to be regular if the Hessian matrix 

(wab QqAQqB ^ 

is regular, and in this case, ul is a symplectic form on TQ. 

We define the Legendre transformation as a fibred mapping FL : 
TQ — y T*Q such that 

(FL(v q ),0 q ) = {X Vq ,a L (v q )) 

where Ttq(X v<i ) — v q G T q Q and {3 q G T*Q. In local coordinates we 
get 

FL(q A jq A ) = (q A ,pA), 
and L is regular if and only if FL is a local diffeomorphism. 

If L is regular, then there exist a unique vector field £l on TQ sat- 
isfying the symplectic equation 

i (L u L = dE L (3.1) 

and moreover, it automatically satisfies the second order differential 
equation (SODE) condition, i.e. 

S£ L = A. 

If, in addition, we assume that L is hyperregular, that is, FL : 
TQ —y T*Q is a global diffeomorphism; then we can define a (global) 
hamiltonian function h : T*Q — y ~R by h = El ° FL' 1 . It is easy to 
show that FL*ujq = ujl and that £l and, then Xh are FL-related. So, 
the solutions of the Euler-Lagrange equations transform by FL into 
solutions of the Hamilton equations and viceversa. Given a vector field 
Z on Q we define a new vector field on Q by 

tl = Tt q o£ l oZ , 
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that is, we have the following commutative diagram 



TQ 



T(TQ) 



Ttq 



Q 



el 



TQ 



Now, Theorem 12.21 can be reformulated as follows. 

Theorem 3.1. Let Z be a vector field on Q such that FL o Z is a 
closed 1-form. Then the following conditions are equivalent: 

(i) £f and £l are Z -related; 

(ii) d(E L oZ) = 

Proof: The result follows as a direct consequence of Theorem 12.21 
and the fact that £l and Xu are FL-related. □ 



4. The Hamilton- Jacobi theory for singular lagrangian 

SYSTEMS 

In this section we shall give a geometric approach to the Hamilton- 
Jacobi theory in terms of the Gotay-Nester-Hinds constraint algorithm 

HOT- 

Let L : TQ — > R be a singular lagrangian, that is, the Hessian 
matrix 

\Wab QqAQgB ^ 

is not regular, or, equivalently, the closed 2-form ujl is not symplectic. 
Therefore, the equation 

it=oj L = dE L (4.1) 

has no solution in general, or the solutions are not defined everywhere. 
Moreover the solutions do not necessarily satisfy the SODE condition. 
Recall that SODE condition is 

S£ = A (4.2) 

or, equivalently, 

Tt q (X) = T TQ (X), 

where tq : TQ — > Q and t T q : TTQ — > TQ are the canonical projec- 
tions. 

Singular lagrangian system have been extensively studied by P.A.M. 
Dirac and P. Bergmann (seedirac), in order to obtain a procedure for 
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canonical quantization of local gauge theories. They developed an al- 
gorithm (called Dirac-Bergmann theory of constraints) that has been 
later geometrized by M.J. Gotay and J. Nester [HI [T2] . 

In the sequel, we will assume that L is almost regular, which means 
that: 

• Mi = FL(TQ) is a submanifold of T*Q; 

• The restriction of the Legendre mapping FL\ : TQ — > Mi is 
a submersion with connected fibers. 

In this case, Mi is called the submanifold of primary constraints. 

If L is almost regular, since ker(TFL) = ker(co>L) D V(TQ), where 
V(TQ) denotes the vertical bundle, and the fibers are connected then 
a direct computation shows that projects onto a function 

hi : Mi — y R . 
Denote by ji : Mi — > T*Q the natural inclusion and define 

Consider now the equation 

i x oji = dhi . (4.3) 

There are two possibilities: 

• There is a solution X defined at all the points of Mi; such X is 
called a global dynamics and it is a solution (modulo ker u>i). 
In other words, there are only primary constraints. 

• Otherwise, we select the submanifold M2 formed by those points 
of Mi where a solution exists. But such a solution X is not 
necessarily tangent to M2, so we have to impose an additional 
tangency condition, and we obtain a new submanifold M 3 along 
which there exists a solution. Continuing this process, we ob- 
tain a sequence of submanifolds 

• • • M fe • • • ^ M 2 ^ Mi ^ T*Q 

where the general description of M; + i is 

M + i := {p G Mi such that there exists X v G T, p M\ satisfying ix^i = dhi}. 

If the algorithm stabilizes at some k, say M^+i = .\ //,.. then 
we say that M^ is the final constraint submanifold which is 
denoted by Mf, and then there exists a well-defined solution 
X of (USD along M f . 

Remark 4.1. There is another characterization of the submanifolds 
Mi that we will describe now. If iV is a submanifold of Mi then we 
define 

TN L = {Z G T p (Mi), peN such that u x (X, Z) = for all X G T p N}. 
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Then, at any point p G Mj there exists X p e T p M; verifying ix^i = 
dh x if and only if (TM fc \ dh x ) = 0, (see [To]). 

Hence, we can define the / + 1 step of the constraint algorithm as 
M t+1 := {p E Mi such that (TM^,dhi)(p) = 0} , 
where TM\~ is defined as above. o 

4.1. Case I: There is a global dynamics. In this case there exists 
a vector field X on M\ such that 

(zjcWi = dhi)\ Ml . 

Moreover, we have tti(Mi) = Q, where tti is the restriction to M\ of 
the canonical projection tiq : T*Q — > Q. 

Next, assume that 7 is a closed 1-form on Q such that j(Q) C Mi. 
Define now a vector field X 7 on Q by putting 

X 7 = T717 0X07. 

The following diagram summarizes the above construction: 

TMi 



T7V1 




Here 71 denotes the restriction of 7. 
We have 

T*(*x-r 71 (X7) wi) = 7i*(«x wi) - 7*(*r 7l (X7) Wi) 

= 7l^(^l) - 7*(«T 7 i(X7) Wi) 

= d^oTi) 
since 7* (ir^jn 07) = ijsn (7*^1) = 0, because 

7*^i = IiJiWq = (ji o 7i)*w Q = 7*w Q = -<fry = 

Therefore, taking into account that Vtti © Tji(TQ) = TM X and loi 
(as it happens with ojq) vanishes acting on two vertical tangent vectors 
with respect to the canonical projection 717 : M x — > Q, we deduce the 
following: 

X - T 7l (X 7 ) e ker(wi) o 7l ) = (4.4) 

Moreover, we will show that it is possible to refine condition ( 14. 4 p 
and to prove that X and X 7 are 71 related. 

First at all, it is clear that for any point p of Mi 
T p (T*Q) = T p Mi + V P (T*Q) 
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where V(T*Q) denotes the space of vertical tangent vectors at p. 

In addition, X — T 7 i(X 7 ) is vertical at the points of Im(7i), so given 
any Z G V P (T*Q), p G Im(7i), we deduce 

u Q (X - T 7l (X 7 ), Z) = along Im( 7l ) 

since ujq vanishes on two vertical tangent vectors. 
Now, given Z G T V M\ we have 

u Q (X - t 7i (x 7 ), Z) = Ul (X- T 7l (X 7 ), Z) = 

because X—T 7 i(X 7 ) G ker(wi), and we obtain that wq(X— T 7i (X 7 ), Z) 
= for any tangent vector Z G T p (T*Q) on any point of Im( 7l ). Since 
ojq is non-degenerate we deduce that X = T 7l (X 7 ) on Im( 7l ). 

In conclusion, we have the following result 
Proposition 4.2. 

X and T 7l (X 7 ) are 7l related <^ d{h\ o 7l ) = 

Remark 4.3. As a consequence of the above result, if hi is constant 
along 7 i(Q) then 7l maps the integral curves of X 7 on integral curves 
of X. So o 7l ) = can be considered as the Hamilton- Jacobi 
equation in this case. o 

4.2. Case II: There are secondary constraints. In this case, the 
algorithm produces a sequence of submanifolds as follows 

• • • M k H • ■ ■ H M 2 ^ Mi ^ T*Q 

We assume that the projections Q r := iiQ(M r ) are submanifolds, and 
that the corresponding projections 7r r : M r — > Q r are fibrations, where 
7r r is the restriction of ttq to M r . 

The constraint algorithm produces a solution X of the equation 

(i x wi = dhi)\ Mf , 
where X is a vector field on Mf. 
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Coming back to the Gotay-Nester-Dirac algorithm we can summarize 
the situation in the following diagram: 



TQ T*Q 



FLi 



M l 



M- 



TTl 



Q 



2 > Q2 



M f Q f 

Assume now that 7 is a closed 1-form on Q such that 

• 7 (Q) C Ml 

• 7 (<2/) C Mf. 

As in Case I, 7 allows us to define a vector field X 1 e £((?/) by 

X 7 = T7T/ 0X071 . 



M 



.f 



7/ 



XT 



TQf 



Qf — 

Here 7/ is the restriction of 7 to Q/- 
Now, given q & Qf, we have 

i (X( 7l (g))-T 97/ (X7( g) )) ^1 ° T <?7l = «X(7(9)) ^1 O T ?7 i - i T , 7/ (X7(,)) Wl o T g7 i 

= dhdfiq)) o T 9 7 = o 7l )(g) . 

Observe that since 7/ is the restriction of 71 we have T g 7/(X 7 (g)) = 
T, 7 (X 7 (g)). Therefore, given F 9 G T,Q then T g7l (F(g)) = T ql {Y{q)), 
and we deduce that 

l T ?7/ (^( 9 ))Wior 57 (y( g )) = Ul (T, 7 (P( 9 )) ) T 97 (y( g ))) 

= (7 1 ^ 1 )(X 7 (g),F(g)) = 7l J> Q (X 7 (g), Y(q)) = d 1 (X^(q),Y(q)) = . 

The previous discussion can be applied to every point q G Qf, therefore, 
taking into account that ou 1 vanishes acting on two vertical tangent 
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vectors and Vn 1 © Tji(TQ) = TMi, we can deduce the following 
X - T 7/ (X 7 ) e ker(wi) & d{h x o 7l )| Q = 

M f 

Usin g a similar argument that in Case I, it is possible to deduce 
that X and X 1 are 7/ related since we have 

T P (T*Q) = T V M X + V P (T*Q) 

for all p e M f . 

Therefore, we deduce the following. 

Proposition 4.4. 

X and T7i(X 7 ) are 7/ related o 7 i)|q / = 

Remark 4.5. Notice that the condition X — T 7 j(X 7 ) 6 Ker(coy) along 
lf(Qf) implies that, if a : R — )■ is an integral curve of X 7 , then 
<7 7 = 7 o cr : R — y Mf is an integral curve of X 

Therefore, the condition 

d(h 071) !Q/ = 

could be still considered as the Hamilton- Jacobi equation in this con- 
text, o 

4.3. Hamilton-Jacobi theory for further geometric constraint 
equations. Besides of the equation ixOJ = dh\ introduced in the pre- 
vious section, other equations have been studied in the literature start- 
ing with the same data, that is, a singular lagrangian. For complete- 
ness, in this section we will discuss some of these equation of motions. 

A good reference for these topics is [9]. 

4.3.1. Extended equation of Motion and the Dirac conjeture. A con- 
straint is called first class provided its Poisson bracket with every other 
constraint weakly vanishes, and second class otherwise (see Section 14.41 
for more details). Dirac [8 J conjectured that all first-class secondary 
constraints generate 'gauge transformations' which leave the physical 
state invariant. See, for instance [10J and references therein for the 
discussion about the avaibility of the Dirac conjeture. Moreover, the 
motivation of our study will be more clear in section 14.41 

Withour entering in physical discussions, we will analyze if it is possi- 
ble to extend our Hamilton-Jacobi formalism for the equations derived 
assuming Dirac conjeture. Therefore, we need first to discuss the ge- 
ometry of this 'extended equation' for singular lagrangians. 

Suppose that we are in the conditions of the previous section. We 
have ji : M\ — > T*Q where M\ is a submanifold and ji the inclusion, 
and a constrained hamiltonian hi : M\ — > R. As before, we study the 
presymplectic system (M 1 ,ui,dhi) and apply the Gotay-Nester-Hinds 
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algorithm, assuming that we reach to a final constraint manifold Mf. 
Denote by jf : Mf — > Mi the inclusion. Now we say that a vector field 

X on Mf is a solution of the extended equations of motion if X can be 
writen 

X = Y + Z (4.5) 

where Y and Z are vector fields on Mf, such that iybj\ = dh\ and 
Z G ker(w/) where ujf := jj(u)x). 

We can now obtain a less restrictive version of the previous Hamilton- 
Jacobi theory, which gives solutions of the extended equations of mo- 
tion. 

Assume again that 7 is a 1-form on Q such that 

(i) Im( 7 ) C Mi 

(ii) Im(7/) C Mf 

(iii) dj = 

From a fixed solution X of the extended equation, we can define 

X I = T-Kf 0X0 jf. 

Proceeding as in the previous section, we have 

l*f{ix-Tif{X"t) w/) = l}{ix w/) - l}{i Tlf {x^) w/) 

= d{h x ojf). 

since 7 is closed. 

Using similar arguments that in the previous section, we deduce the 
following 

Proposition 4.6. Under the above conditions, we have 
X - Tjf(X^) E ker(w / ) <S=> o 7/ ) = 0. 

M f 

Ifd(h x o 7/ ) = 0, then X - Tj f (X^) e kex(u f ). 

Proof: It follows the same lines of the proofs of Proposition 14.41 but 
now observing that 

TM f = Tjf(TQ f ) © V7T f 
and Viif C V M; ^q- □• 
By the last proposition Tjf{X 1 ) = X + Z, with Z G kerujf. Then, 
form (g3D we have that Tj f {X^) = Y + (Z + Z). So, Tj f (X^) is a 
solution of the extended equations of motion. 

Therefore, the condition 

d{h\ o jf) = 

could be still considered as the Hamilton- Jacobi equation in this con- 
text. 
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4.3.2. Hinds algorithm. Besides of the Gotay-Nester-Hinds algorithm, 
other approaches have been discussed in the literature. In particular 
we briefly recall the algorithm introduced by Hinds (see Gotay [U] 
for a detailed discussion). Hinds algorithm also start considering the 
equation ix^i = dhi as the Gotay-Nester-Hinds algorithm. The algo- 
rithm generates a descending sequence of constraint submanifolds. In 
the favorable case, the algorithm stabilizes at a final constraint sub- 
manifold which we will denote again by Nf (see discussion below). It 
is important to point out that, in general, this constraint sub manifold 
Nf will be different from the final constraint submanifold obtained by 
the Gotay-Nester-Hinds algorithm, that is Nf ^ Mf. In principle, both 
algorithms start to diverge from each other after the second step. 

In more geometric terms, assume that we are in the conditions of the 
previous section. Define Ni := Mi as we did before and denotes A^ +1 
for I > 1 the following subset 

iV; +1 := {p G Ni such that exists X G T p Ni verifying ix^i = dhi}, 

where, if we call ki : Ni — > N± the natural inclusion, then ui := fc z *u;i 
and hi := k\h\. We obtain the sequence of submanifolds 

• • • N k ■ ■ ■ N 2 iVi = M x T*Q. 

Again if the algorithm stabilizes, i.e. = iV&. +1 , then we say that N k 
is the final constraint manifold, Nf. In this case, the Hinds algorithm 
produces a solution X G %(Nf) of the equation 

ix^f = dhf. 

This equation is less restrictive than (14.3)) . and so the two algorithms 
diverge for I > 2. We will come back later to the above equation. 

Now, we can develop a Hamilton- Jacobi theory in this setting. 

Assume that there exists a 1-form 7 on Q satisfying 

(i) Im( 7/ ) C N f , 

(ii) g?7 = along Nf. 

Then we can define X 1 = Titf o X o 7^ and state the equivalent 
Hamilton- Jacobi theory. The proof follws the same lines that in propo- 
sition 14.61 

Proposition 4.7. 

X - T 7/ (X 7 ) G ker(w/) d(h o 7/ ) = 0. 

4.4. Relation to the Dirac-Bergmann theory of constraints. In 

this section we will discuss the relation of the Gotay and Nester theory 
with the original Dirac-Bergmann theory of constraints. 

Assume that we begin with an almost regular lagrangian L : TQ — > 
K. Then there exists an open neighbourhood, U C T*Q where in 
canonical coordinates (q A ,PA), Mi D U is given by the vanishing of 
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functions ,Pa) defined on U . The functions $ l are called primary 
constraints. 

Remember that we can project El to h\ : M\ — > 1R, and any exten- 
sion of hi to U should be of the form 

H = h + Ui & , 

where h is an arbitrary extension to U of hi. The functions Ui, 1 < % < 
2 dim Q — dim Mi are Lagrange multipliers to be determined. 

Acording to Dirac the equations of motion are 

•A _ 8H_ _i_ 

V dp A 1 dp A 

8H _ d& 

PA - dqA u i dq A 

which must hold over U\ := Mi D U. If we denote ji : U\ — > U the 
inclusion, and u\ = j*U)q, the preceding equations can be equivalen 
rewritten as 

ix^i = dhi or {ix^Q = dh + UidQ 1 )^ (4.6) 

which are the equations that we have considered in the Gotay-Nester- 
Hinds algorithm. 

Since X must be tangent to U\ we should have 

= (X(^))| C / 1 = {^,^}| 1 / 1 = {^^ + ^^'}| 1 / 1 
= {{¥,h} + Uj {¥,&}) m . 

These equations can be trivially satified, determine some Lagrange 
multipliers or add new constraints on the variables q A , pa over U\. 
These new constraints, if any, are called secondary constraints. Sup- 
pose that we have obtained the secondary constraints £ a . So, we have 
to restrict the dynamics to Ui '■— U\ PI (^ a ) _1 {0}- 

Again, the solution must be tangent to Ui and it requires that 

= (X(^)) ]U2 = {^,H} ]U2 = {^,h + u t ¥} ]U2 

= ({^ /»} + «*{£", * 4 }V a . 

As before, these equations may determine more Lagrange multipliers or 
add new constraints to the picture, that is, new secondary constraints. 
Iterating this procedure, if the algorithm stabilizes, we arrive to a set 
U f which is an open subset of the final constraint manifold Mf obtained 
by the Gotay-Nester-Hinds algorithm (see [T5] for a proof ). 

It is necesary to introduce some definitions. We say that a function 
defined on U is first class if its Poisson bracket with every constraint 
(primary and secondary) vanishes. Otherwise, it is said to be of second 
class. 
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We can reorder constraints into first class or second class. We will 
denote by \ a an d £, b , the primary first and second class constraints, re- 
spectively; and by ip c and 9 d , the secondary first and second class con- 
straints, respectively We will also denote by // , A& the corresponding 
Lagrange multipliers for the primary first and second class constraints, 
respectively. 

So, if the problem has a solution, we must obtain a vector field X 
over Uf, which satisfies the equations 

(i x ujq = dh + fi a dx a + hd£ b )\M r 

The Ab's are determined functions and the // a 's can be varied to obtain 
other admissible solutions. In consequence, it is also clear that primary 
first class constraints correspond to gauge transformations which leave 
the physical state invariant. As we have discussed before, Dirac con- 
jectured that the first class secondary constraints may also generate 
gauge transformations, therefore, the generalized equations of motion 
discussed in Subsection 14.3.11 are locally rewritten as 

(i x u Q = dh + fi a d X a + \ b d? + v c dijj c )\ Mr (4.7) 

where A& are still determined functions and fi a and v c can be varied ar- 
bitrarily. The hamiltonian h + fi a x a + \b£, b + v c tp c is called the extended 
hamiltonian, and equation (14. 7p the extended equation of motion fol- 
lowing the notation of [9J. Geometrically the solutions of (14. 7p are 
just 

X = Y + Z, 

where Y is a vector field on Mj solution of the equations of motion, 
14.61 an d Z e ker(ojf) where Uf is the restriction of uj\ to Mf. 

Remark 4.8. If we proceed in the same way with the Hinds algorithm 
developed in Subsection I4.3.2[ we will arrive to solutions X satisfying 

(i x ujq = dh + /J> a dx a + Abd£ 6 + Vcip c + w-^ d )\ Np 
where t%, w-j are the Lagrange multipliers corresponding to the con- 
straints ijj c and Q d . Note that iff and (f now correspond to the sec- 
ondary constraints of the final constraint manifold Nf in the Hinds 
algorithm. o 

4.5. Examples. Now we illustrate the previous propositions with sev- 
eral examples. 

4.5.1. There are only primary constraints. 

Example 4.9. This example is discussed by O. Krupkova in [16j. Let 
L be the Lagrangian L : TR 3 — > R given by 

L( q \ q \ q \q\e^) = \(e+?y. 



Ker(o>i) 
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Then FL is given by FL : TR 3 ->■ T*R 3 

FL(q\ q 2 , q\ q\ q\ q 3 ) = (q\ q 2 , q 3 , q l + q\ q 1 + q\ 0), 
and the primary constraints are 

^ 1 (q A ,Pa) = Pi ~ P2 $ 2 (q A ,p A ) = Ps- 

So 

Mi = {(q\q 2 ,q 3 ,Pi,P2,P3) € M 6 such that p x = p 2 , p 3 = 0} 
and we can use (q 1 , q 2 , q 3 , Pi) as coordinates on M±. 
It follows that 

El = (q 1 + W + (q 1 + ?)? — L — L 
h = \{pif 

u>i = dq 1 A dpi + dq 2 A dp\ 
d d d 
dq 3 ' dq 1 dq 2 
and a particular extension of the hamiltonian is 

h(q A ,p A )= 1 -( Pl ) 2 

It is easy to see that at the points of M x := Im(FL) 
{¥, h + u 1 ® 1 + m 2 $ 2 } = 0, i = 1,2. 

So we have global dynamics on Mi and it holds that 

| $ 1 $ 2| = Q 

and we conclude that there are only first class constraints. The solu- 
tions of (ixWi = dhi)\M 1 on Mi are given by 

d d d d 

^ 1 dq 1 ^ dq 3 ^ 2 dq 1 dq 2 ^ 
where fi and / 2 are functions on M x . 
We now look for 7 e A 1 (Q) such that 

(i) 7 (Q) C Mx 

(ii) o 7) = 

(iii) ^7 = 

Suppose 

7(?\ g 3 ) = (?\ q 2 , q 3 , 7i(?\ q 2 , q')Mq\ q\ q")Mq\ q\ q 3 )), 
then, 7(<5) C Mi implies that 71 = 72 and 73 = 0. 

The condition d(hi 07) = implies that |(7i) 2 = constant, and 
because of that 71 = c where c is a constant and so 72 = c. 

Now 7(g 1 , q 2 , q 3 ) = (q 1 , q 2 , q 3 , c, c, 0) and d^y = is trivially satisfied. 
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If we take the general solution Pi-^t + fi^s + .M^r — ^2) then we 
obtain X y = Tn q oXo 1 = c^ T + (f 1 o 7 )£ + (f 2 7 )(^ - If 

we apply T r y(q A )(X' 1 (q A )) = X( 7 (g A )), then we recover the solution X 
over the points of 7. It is clear, that integral curves of X 1 are applied 
by 7 into integral curves of X along Im 7. 

Example 4.10. This example has been discussed by J. Barcelos-Neto 
and N.R.F. Braga [3J. Let L be the Lagrangian L : TM 4 — > K given by 

L(g\ g 2 , g 3 , g 4 , q\ q\ q\ g 4 ) = tf+f)<£+£?+l (tf? ~ 2?V - (g 3 ) 2 ) . 

Then FL is given by FL : TM 4 ->■ T*M 4 

FL(q 1 , q\ q 3 , q\ q\ q\ q 3 , g 4 ) = (g\ q 2 , q 3 , q\ q 2 + q 3 , 0, g 4 , 0) 
and the primary constraints are 

® l {q A ,PA) =Pi~q 2 - q 3 , $ 2 (q A ,p A ) = P2, 
<f> 3 (q A ,p A ) =Pa- q 4 , <5> A {q A ,p A ) = Pa- 

So 

Mi = {(q 1 ,q 2 ,q 3 ,q A ,pi,p2,P3,P4> £ such that 

Pi = q 2 + q 3 , P2 = 0, p 3 = g 4 , p 4 = 0}. 

and we can use (q 1 , q 2 ,q 3 , g 4 ) as coordinates on M\. 
It follows that 

E L = (q 2 + q^q 1 + g 4 g 3 — L = -±((g 4 ) 2 - 2g 2 g 3 - (g 3 ) 2 ) 

^i = -|((g 4 ) 2 -2g 2 g 3 -(g 3 ) 2 ) 

Ui = dq 1 A dq 2 + dq 1 A <ig 3 + dq 3 A rfg 4 

Ker(wi) = {0} 

so (Mi,ui) is a symplectic manifold. 

It we prefer to follow the Dirac-Bergmann algorithm, then one should 
take an extension h(q A ,pA) = — |((g 4 ) 2 — 2g 2 g 3 — (g 3 ) 2 ) of hi. It is easy 
to see that at the points of Mi := Im(FL) 

h + u 1 ® 1 + u 2 $ 2 + m 3 $ 3 + m 4 $ 4 } = — u 2 - u 3 
{$ 2 , h + u 1 ® 1 + u 2 $ 2 + m 3 $ 3 + m 4 $ 4 } = -g 3 + u l 
{$ 3 , h + u 1 ® 1 + m 2 $ 2 + m 3 $ 3 + m 4 $ 4 } = -g 2 - g 3 + u 1 - w 4 

{$ 4 , /l + U 1 ® 1 + M 2 $ 2 + M 3 $ 3 + M 4 $ 4 } = It 3 + It 4 , 

which determine completely the Lagrange multipliers: 
u 1 = g 3 , u 2 = g 4 , u 3 = -g 4 , u A = -g 2 , 
and then all the constraints are of second class. 

The solution of the equation (ix^i = dhi)\ Ml is given by 

o 8 A 8 A 8 , n n. n 
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We will study now the solutions of the Hamilton- Jacobi equation. 
So, we look for 7 e A 1 (M 4 ) such that 

(i) 7 (R 4 ) C Afi 

(ii) d(hx o 7 ) = 

(iii) d 7 = 

If 7 (g A ) = (?' 4 , 7i(g A ), 72(^)5 73(g^),74(g A )), then the condition 
j{Q) C M l gives 

7i(g 1 ,g 2 ,g 3 ,,g 4 ) = g 2 + g 3 
7 2 (g 1 ,g 2 ,g 3 ,,g 4 ) = 

73<y,g 2 ,g 3 , ,g 4 ) = g 4 
7 4 (g 1 ,g 2 ,g 3 ,,g 4 ) = 

But ^07= — |((g 4 ) 2 — 2g 2 g 3 — (g 3 ) 2 ), so the equation d(h 1 07) = 
if and only if 7(g) = {q A , 0, 0, 0, 0). 

Example 4.11. This example has been discussed by K. Sundermeyer 
[30J. Let L be the Lagrangian L : TR 2 — > R given by 

L( q \ q ^q\f)= 1 -( q y+eq 1 +q 2 <l 1 . 

Then FL is given by FL : TR 4 -»■ T*R 4 

FL(q\q 2 ,q\q 2 ) = (q\q 2 ,q l + q 2 ,q 1 ) 
and the primary constraints are 

®\(1 A ,Pa) =P2-q 1 

So 

Mi = {(g 1 , q 2 ,Pi,P2) G R 4 such that p 2 = g 1 }, 
and we can use (gSg^Pi) as coordinates on M\. 
It follows that 

El = \q X 

hi = \{pi~q 2 ) 

Ui = dq 1 A dp 1 + dq 2 A dq l 

KerM = (wr w) 

Let i 

h(q A ,PA) = ^(pi- g 2 ) 

be an extension of the hamiltonian. 

It is easy to see that at the points of Mi := lm(FL) 

+ = 
and therefore we have global dynamics. 
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The solution of the equation (ix^i = dhi)\Mi ls gi ven by 

x = (pi - ? 2 ) 4i +f4-2+f^-+ (pi - ? 2 )^> 

ag 1 ag^ api ap 2 

where / e C°°(Mi) 

If we now look for 7 G A 1 (M 4 ) such that 

(i) 7 (M 4 ) C Mi 

(ii) d{h x o 7) = 

(iii) c/7 = 

then 7(g 1 , g 2 ) = (g 1 , g 2 , 7 1 (g 1 , g 2 ), 7 2 (g 1 , g 2 )) given by 

l{q\q 2 ) = {q\q\q\q 1 ) 

satisfies all the requiered conditions, because Pi(7(g 1 ,g 2 ) = g 2 , 7 = 
d^q 1 • g 2 ) and hi o 7(g 1 ,g 2 ) = \{q 2 — g 2 ) = 0. Given an arbitrary 
solution X = (p 1 -q 2 )-^ T + f-^ + f-^- + (p 1 -q 2 )-^-oi the constarined 
dynamics, we have that 



•V" (/°7)v, 



d_ 
dq 

and also 

r 7 (jro = (/o 7 )A + (/07 )| r 

which is precisely X along Im(7). 

4.5.2. There are secondary constraints. 

Next, we are going to describe several examples where secondary 
constraints appear. 

Example 4.12. This example has been discussed by M.J. Gotay and 
J.M. Nester [14] . Let L be the Lagrangian L : TR 2 ->■ R given by 

L(q\q 2 ,q\q 2 )= 1 -(q 1 ) 2 + q 2 (q 1 ) 2 - 
Then FL is given by FL : TR 2 -»■ T*M 2 

FL(g 1 ,g 2 ,g 1 ,g 2 ) = (g 1 ,g 2 ,g 1 ,0) 
and the primary constraints are 

$ x (g A ,p A ) =P2- 

So 

Afi = {(q\q 2 ,Pi,P2) e M 4 such that p 2 = 0} 
and we can use (gSg 2 ,^!) as coordinates on M x . 
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It follows that 

^ = |(g 1 ) 2 -g 2 (g 1 ) 2 

hi = \{p 1 ?-q 2 {q 1 ) 2 
ui\ = dq 1 A dpi 

KerM = (£) 

Let 

h(q A ,p A ) = \(pi) 2 -q 2 (q 1 ) 2 

be an arbitrary extension of the constarined hamiltonian hi to T*IR 2 . 
It is easy to see that at the points of M x := Im(FL) we have 

+ = -(g 1 ) 2 

and therefore we need to restrict the dynamics adding a new constraint 

&(q A ,PA) = q 1 . 

Now M 2 := {(q 1 ,q 2 ,Pi,P2) e K 4 such that p 2 = 0, q 1 = 0} and Q 2 := 
tt q (M 2 ) = {(q l ,q 2 ) e M 2 such that q 1 = 0}. We have on M 2 

{Q^h + u 1 ® 1 } = 
{^ 2 ,h + u 1 ^ 1 }=p 1 , 

and we need to restrict again the dynamics, adding the constraint 

$ 3 (? A ,Pa) =Pi 

Now M 3 := {(q 1 ,q 2 ,pi,p 2 ) e such that p 2 = 0, g 1 = 0, p 1 = 0} 
and Qs = Q 2 = {(q 1 , q 2 ) G M 2 such that q 1 = 0}. Along M 3 we have 

+ = 

{$ 2 ,/i + w 1 $ 1 } = 
{^/i + u 1 ^ 1 } = 0, 

and M3 is the final constraint sub manifold, Mf. We can easily check 
that "I) 1 is a first class constraint and $ 2 , $ 3 are second class. 

The solutions of the equation (ix^>i = dhi)\M 3 are of the form 

X-f 9 

where / G C°°(M 3 ). 

A solution of the Hamilton- Jacobi equation, should be 7(? 1 ,? 2 ) = 
(g 1 ,g 2 ,7i(g 1 ,g 2 ),7 2 (g 1 ,g 2 )), such that 

(i) 7 (Q) C Mi and 7/ (Q / ) C M, 

(ii) d(h\ o 7)|Q/ = 

(iii) ^7 = 
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The condition j(Q) C Mi implies 72 = 0. Next we compute d'y, 

dl = ^l dq 2 A dgi + ^g 1 A dq 2 = ^ dq 2 a dq l = 
c>g 2 ag 1 ag 2 

and we deduce that must vanish and 71 must be a function of g 1 . 
The condition d(/ii 07)^ = can also be easily computed. We have 

d(hi o 7) = c /(I( 7l ) 2 - gV) 2 ) = (71 1^ - 2gV)dg 1 + (g 1 ) 2 ^ 2 

and, along Qf = {(g 1 , g 2 ) G M 2 such that g 1 = 0}, we deduce 

^io 7 )|Q / =7i(0)|^(0)dg 1 

For example, if we take 71 = g 1 , 72 = 0, all the above conditions are 
satisfied, and 7/(<5/) C Mf. 

Now, take a solution X = f-£p] at the points of Qf we get 



and so 



X^(0,g 2 ) = (*■,).(/((), g 2 , 0,0) A) = /(0,g 2 ,0,0) A 
T 7/ (X^(0,g 2 )) = /(0,g 2 ,0,0)A ) 



and we obtain the solution X along Im(7/). 

Example 4.13. This example has been discussed by M.J. Gotay |10j . 
Let Q := {(g^g 2 ) € M 2 such that g 1 7^ 0} and L be the Lagrangian 
L : TQ — > K. given by 

L(g\g 2 ,g\g 2 ) = A(«p)* 

Then FL is given by FL :TQ ^T*Q 

FL(g 1 ,g 2 ,g 1 ,g 2 ) = (g 1 ,g 2 ,0,g 2 /g 1 ) 
and the primary constraints are 

$\q A ,p A ) = Pl . 

So 

M — {{q 1 iQ 2 £ suc h that pi = 0}. 
and we can use (q 1 ,q 2 ,P2) as coordinates on M x . 
It follows that 

E L = L 

hi(q A , PA ) = i(P2) 2 
LOi = dq 2 A dp 2 

Ker(wi) = 
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Let h(q A } pA) = \{P2) 2 be an extension of the hamiltonian. 
It is easy to see that at the points of M\ := Im(FL) we get 

and therefore we need to restrict the dynamics adding a new constraint 

® 2 (q A ,PA) =P2- 

Now M 2 := {(q L ,q 2 ,Pi,P2) £ TQ such that p± = 0, P2 = 0} and Q2 '■— 
7Tq(M2) = Q. At the points of M2 we have 

+ = 

{Q^h + u 1 ® 1 } = 0, 

and M 2 is the final contraint manifold. From $ 2 } = we deduce 
that the constraints are all first class. 

The solutions are of the form X = /^r where / G C°°(M2). 

If we look for a solution of our Hamilton- Jacobi equation, 7, such that 
7(g\g 2 ) = (g 1 ,g 2 ,7i(g 1 ,<? 2 ),72(<? 1 ,<? 2 )), then the condition ^(Q) C M 1 
implies 7 = and 7/ = 0. All conditions are verified and, given a 
solution X, we obtain that X 1 and X are trivially 7^-related. 

Example 4.14. This example has been discussed by R. Skinner and 
R. Rusk [29J. Let L be the Lagrangian L : TM 3 — > K given by 

L(g\ g 2 , g 3 , g\ g 2 , g 3 ) = ^g 2 (g 3 ) 2 + g'g 3 . 

Then FL : TM 3 T*R 3 is given by 

FL{q\ g 2 , g 3 , g 1 , g 2 , g 3 ) = (g 1 , g 2 , g 3 , g 3 , 0, g 1 ) , 

so that we have a primary constraint Q l (q A 1 pA) — Pi- This means that 
the primary constraint submanifold is 

Mi = {(q\q 2 ,q\pi,P2,p 3 ) G R 6 such that p 2 = 0}, 
and then we can use (g 1 , g 2 , q 3 ,Pi,P3) as coordinates on M x . 
It follows that 

El = g 3 g* + g*g 3 — L = — |g 2 (g 3 ) 2 + g : g 3 

hi(q A ,p A ) =pm - y 2 (q 3 ) 2 

loi = dq 1 A dpi + dq 3 A dp3 

Ker( Wl ) = (^)- 

As in the previous cases, take an arbitrary extension of the hamiltonian 
hi, for instance 

h(q A ,PA) = pm - 2<? 2 (<? 3 ) 2 - 

It is easy to see that at the points of Mi := lm(FL) 

+ = |(g 3 ) 2 
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and therefore we should restrict the dynamics adding a secondary con- 
straint 

® 2 (q A ,PA) = q 3 - 

Now M 2 := {(q 1 ,q 2 ,q 3 ,Pi,P2,P3) G ^ 6 such that p 2 = 0, g 3 = 0} and 
Q 2 ■= ttq(M 2 ) = {(g 1 ,? 2 ,? 3 ) G M 3 such that g 3 = 0}. Along M 2 , we 
have 

+ = 

{$ 2 , /i + m 1 ^ 1 } = p 1 . 

Therefore, we need again to restrict the dynamics, adding the con- 
straint <& 3 (q A ,pA) = Pi .NowM 3 := {(g 1 , q 2 , g 3 ,Pi,P2,£>3) £ ^H 6 such that 
0, g 3 = 0, Pl = 0} and Q 3 = Q 2 = {(g 1 ,^ 3 ) G K 3 
such that g 3 = 0}. Along M 3 we have 

+ = 

{S^/i + w 1 * 1 } = 
{$ 3 ,/i + w 1 $ 1 } = 

and then M 3 is the final contraint manifold, denoted by My; therefore, 
Qf = Qs- We deduce that the constraints are all first class. 
The solutions of the equation (ix^i = dhi)\ M , A are of the form 

x-p 8 + / 8 

dq 1 dq 2 ' 

where / E C°°(M 3 ). 

Now we look for a solution of our Hamilton- Jacobi equation, that is 
7(g 1 ,g 2 ,g 3 ) = (g 1 , g 2 , g 3 , 7i(g 1 , g 2 , g 3 ), 7 2 (g 1 , g 2 , g 3 ), -f 3 (q\ q 2 , g 3 )), such 
that 

(i) 7 (g) C Mi and 7/ (Q / ) C M f 

(ii) d{h 1 o 1 ) ]Qf = 

(iii) d-f = 

The condition 7(Q) C Mi implies 7 2 = 0; the condition 7/(Q/) C My 
implies (7/)i = for i = 1, 2 and, the condition d(/ii o 7) = is 

#i°7) = iiis - \q 2 {q 3 ) 2 

= (gHs + f?7i) dg 1 + (|J73 + f?7i + dq 2 

+ (&7a + fpbi + ?V) ^ 3 

Hence, 

7)|Q, = (&73 + &7i) <V + (&73 + g§7i) ^ 2 
+ (&73 + gl7i) dq 3 



26 M. DE LEON, J. C. MARRERO, D. MARTIN DE DIEGO, AND M. VAQUERO 

The condition g?7 = implies 

dl= ^dq 2 Adq 1 + ^dq 3 Adq 1 

+ ^dq 1 A dq 3 + fjjdg 2 A dq 3 = 

taking into account that 72 = 0, and therefore 

971 _ n 
a 9 2 ~ u 

ch/i _ 873 
dq 3 dq 1 

973 _ n 

a 9 2 ~ U 

A particular solution is obtained putting 71 = 72 = 0, and 73 an 
arbitrary function of g 3 , for example 73 = q 3 . 

For instance, take X = P3^t + and S = |(g 3 ) 2 , then 

1 (q\q 2 ,q 3 ) = (q\q 2 ,q 3 ,0,0,q 3 ) 
and at the points of Qf we obtain 



d d d 

so that 



X\q\ q\ 0) = T7T f (0—+f(q 1 , q\ 0, 0, 0, 0)—) = f(q\ q 2 , 0, 0, 0, 0)—, 



T ls {X\q\ q 2 , 0)) = f{q\ q 2 , 0, 0, 0, 0)^ 

We can also apply proposition 14.61 to the latter example and obtain 
solutions of the extended equation. 

For instance, consider 7(g 1 , q 2 , q 3 ) = (g 1 , q 2 , q 3 , 7i(g ), l2(q A ), l3(q A )) 
given by 

i(q\q 2 ,q 3 ) = (q\q 2 ,q 3 ,q 3 ,o,q 1 ) 

We have 

(i) 7/(<5/) C Mj because 

j f (Q f ) = {(q\q 2 , 0,0,0^^ G R 6 such that q 1 , q 2 E R}. 

(ii) If we take coordinates (g 1 , q 2 ) in Qf, then d{hi 07^-) = d(0 ■ — 
g 2 -0) = 0. 

(iii) c?7 = 0, in fact, 7 = d(q l ■ q 3 ). 

If we consider a solution X = ps^r + / A , we can compute 

X V, g 2 , 0) = Ttq 0X0 7/ (g\ g 2 , 0) = g 1 ^ + f{q\ g 2 , 0)^ 
and also 

T 7/ (*V, ^ 0)) = g 1 ^ + /(g\ g 2 , 0) A + ^ 
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which is a solution of the equation %x w 3 = dh 3 where, if i 3 : M 3 — > T*Q 
is the inclusion on T*Q and j 3 : M 3 — > M 1; then co 3 = ^(wq) and 
Ml':)- 

Note that 7 in this case is not a solution of our Hamilton- Jacobi 
problem because d{h o 7)(g 1 , q 2 , 0) = q 1 dq 3 7^ 

4.6. Relation to classical Hamilton-Jacobi Theory. In this sec- 
tion we will connect the Hamilton-Jacobi theory developed in the pre- 
vious sections with the classical Hamilton Jacobi theory on T*Q using 
an appropriate extended hamiltonian. 

We will use the same notation that in section 14.41 We start with an 
almost regular lagrangian L : TQ — > R, and then lm(FL) = M\ is a 
differentiable submanifold of T*Q and, in addition, we can define hi 
implicitly by hi o FL = El. We denote U\ = j{u)Q, where ji : Mi — >• 
T*Q is the inclusion and uq is the canonical symplectic form of the 
cotangent bundle. We take local coordinates (q A ,p A ) in an open set 
U C T*Q, such that Mi is given locally by the vanishing of independent 
functions & l (q A ,p A ), called primary constraints. 

Remember that the equations of motion have the form (ix^i = 
dhi)\u i: where Ui := MidU. This equations are equivalent to (ix^Q = 
dh + LiidQ 1 )^ where h is any extension of hi to U defined on Mi and 
/ij are Lagrange multipliers. 

4.6.1. Case I: There are only primary constraints. First, we suppose 
that there exist a global solution X, i.e. X is a vector field on Mi that 
satisfies the equations of motion. We reorder constraint functions in 
two classes: first class constraints denoted by \ a an d second class con- 
straints denoted by We also denote by u a and A& the corresponding 
Lagrange multipliers. Then the equations of motion are 

(i x uq = dh + u a dx a + A 6 d£ 6 )|M r 

Now, suppose that u a and Af, are functions defined on U. It is clear that 
X is the restriction to Mi of the hamiltonian vector field corresponding 
to a hamiltonian of the form h + u a x a + K£, b - Ln fact, all the solutions of 
the equations of motion are obtained in this way varying the functions 
u a arbitrarily and with prescribed values of A&. Next, we are looking 
for a solution of our Hamilton-Jacobi problem, that is, a 1-form 7 
satisfying 

(i) d~i = 

(ii) Im(7) C Mi 
(hi) d(h o 7) = 

Condition (iii) can be easily checked that it is equivalent to d((h + 
UaX a + o 7) = because (h + u a x a + A b £ b )| Ml = fa. So, it is 
evident that the solutions of the classical Hamilton-Jacobi equation for 
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the hamiltonians h + u a x a + A b £ fe (where u a are arbitrary functions and 
the rest are fixed) inside M\ and the solutions for our Hamilton- Jacob i 
problem coincide. 

4.6.2. Case II: The general case. Suppose now that the algorithm do 
not stop at Mi, then we obtain the sequence of manifolds 

• • • M k • • • M 2 <-> Mi <-> T*Q. 

and we suppose that the algorithm stabilizes in a manifold Mf of 
dimension> 0. 

We can reorder the constraints in first and second class (maybe 
changing the independent set of constrainsts) . We will denote x a an d 
£ b the primary first and second class constraints and by i/j c and 9 d 
the secondary first and second class constraints. We will also denote 
by u a ,Xb,v c and the corresponding Lagrange multipliers. Again a 
solution X of the equations of motion verifies 

(i x uj q = dh + u a dx a + A b d£ 6 )| M/ . 

As above, X is the restriction to Mf of the hamiltonian vector field 
given by the hamiltonian h + u a x a + Aft£ 6 where some multipliers are 
determined applying the constraint algorithm. 

We are looking for 7 e A X (Q) satisfying 

(i) c/7 = 

(ii) Im( 7 ) C Mi 

(iii) 7(Qy) C M f 

(iv) d(hx o 71) = 

Note that (iv) is equivalent to the equation d((h + u a x a + A&£ 6 ) 07) = 
because (h + u a x a + A^ 6 )^! = hi, and so, the solutions of the 
classical Hamilton- Jacobi theory contained in Mf for the hamiltonians 
h + u a x a + Ab£ b are just the solutions of our Hamilton- Jacobi problem. 

4.7. Relation to other theories. The Hamilton- Jacobi theory for 
degenerate lagrangians have been discussed by several authors in the 
last 20 years. Let us recall some previous attempts. 

(i) In the papers by Longhi et al. |2H [25] it is discussed the 
case of a time independent lagrangian which is homogeneous 
in the velocities. It is shown that we can substitute an ar- 
bitrary lagrangian by an homogeneous one using the tradi- 
tional procedure by adding new variables and, then, this new 
homogeneous lagrangian has zero energy. The authors show 
that the hamiltonian can be added as a new constraint and, 
in consequence, they restrict themselves to the case when the 
hamiltonian is identically zero. The integrability condition for 
the resultant Hamilton- Jacobi equations implies that they can 
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only consider first class constraints. On the other hand, in the 
paper by Rothe and F. G. Scholtz [27J an almost-regular la- 
grangian L(t, q A , q A ) is considered. If the Hessian { -q^k§as )i nas 
rank n — m\ then, the constraint submanifold M\ is locally de- 
scribed by coordinates (q A , p a ), where only , a — mi + 1, . . . , n. 
The remaining momenta p a ; a = 1, . . . , mi are functions of 
t, q A , p a , that is, p a = —f a (t,q A ,p a ) and represent the pri- 
mary constraints (fi a (t,q A ,PA) = Pa + f a (t,<l A ,Pa)- Then they 
consider the system of partial differential equations 

f + M^,i|) = o 

f| + / a (t,g A ,JJ) = 6=1,. ..,mi 

where hi is the hamiltonian defined on the primary constraint 
manifold by the projection of the lagrangian energy. 

(ii) The theory discussed in [5] is similar to our theory in the case of 
global dynamics, but they do not take into account secondary 
constraints. The authors also use the lagrangian homogeneous 
formalism to obtain the standard Hamilton- Jacobi theory for 
time dependent systems. 

(iii) M. Leok and collaborators [UJ use the Dirac structures setting, 
and secondary constraints are not considered. 

4.8. Lagrangian setting. The equations of motion are globally ex- 
pressed by the presymplectic equation 

i i u L = dE L , (4.8) 

where a possible solution £ is not in principle a SODE. 

Therefore, in addition to the problem of finding solutions for (14. 8p . 
we must study the second order differential problem, that is, we shall 
obtain a solution of (I4.8[) satisfying the additional condition St; = A. 

If we apply the constraint algorithm to the presymplectic system 
(TQ,UL,dEi) we obtain a sequence of submanifolds. 

... Pk ^ ...^ P 2 ^ P i:=T Q 

Assume that the algorithm stabilizes at some Pk+i = Pk — Pf, which 
is the final constraint submanifold. 

If we consider, as above, the presymplectic system (Mi, ui, dhi), and 
apply the constraint algorithm to the equation 

%x oji = dhi (4.9) 

we obtain a sequence of submanifolds 

• • • Mfc <->•••'—)■ Mi Mi ■=->■ T*Q. 

such that 

FL(Pi) = Mi, for any i, 
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and 

FU := FL\ Pi iPi^Mt 

are surjective submersions. 

As a consequence, both algorithms stabilizes at the same step, say 
k, and then 

FL(P f ) = M f 

and 

FL f :P } ^M f 

is a surjective submersion. Moreover, we have the following results. 

Proposition 4.15. If I; is a FLf-projectable solution of (14. 8p . then its 
projection T FLf (£) is a solution of (j4.9j) . 

Conversely, if X is a solution of (14. 9p . i/ien an?/ FLf projectable 
vector field on Pf which projects on X, is a solution of (14.81) . 

Next, we shall discuss the SODE problem as it was stated by M. J. 
Gotay and J. Nester [9j [12] (see [I] for an alternative description). 

The results in [H1CI2] can be summarized in the following result. 
Theorem 4.16. 

(i) If!; is a FLf-projectable vector field on Pf then for any p G Mf 
there exists a unique point in each fiber FLj 1 ^), denoted by 
Vt{p) a t which £ is a SODE. The point rj^{p) is given by 

%(?) := Ttq{£(p)) 

(ii) The map 

ft: Mf — ► Pf 

P ->■ Pt(p) '■= Vt(p) 

is a section of FLf : Pf — > Mf and on Im(/3^) there exists a 
unique vector field, denoted by X^, which simultaneously satis- 
fies the equations 

ixt: u L = dE L 
SX^ = A 

We will now recall the construction of a solution of the dynamical 
equation which simultaneously satisfies the SODE condition. If X := 
(FLf)*^), then X is a vector field on Mf satisfying ix uj± = dh\. The 
vector field X^ described in (ii) is given by 

A detailed proof can be seen in [12], but for the sake of complet- 
ness, we recall here the way to choose the points on the fibers as it is 
stated in the Theorem 14.161 (i). 
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In the last part of this section we come back to the Hamilton- Jacobi 
problem, but now in the lagrangian setting. 

The application of the constraint algorithm is summarized in the 
following diagram 




Assume, as before, that Qi = 7Tg(Mj) are submanifolds and 7Tj = 
7TQ| A/ : Mi — > Qi are surjective submersions. Since tq = ttq o FL, 
then tq(Pj) = iiQ^Mi) = Qi, and Pi also projects onto Qi. We denote 

r f = T Q\ Pf ■ Pf -»■ Qf- 

In consequence, the following diagram is commutative. 

Pf 




Now, if X is a solution of ix Wi = dh% on M/ and 7 is a 1-form which 
is a solution of the Hamilton- Jacobi problem, that is, 

(i) y(Q) C Mi and 7/ (Q / ) C My 

(ii) d(/ii o 7l )| Q/ = 

(iii) dy = 

then we can define X 7 = T7170X07J. From Proposition 14.61 we deduce 
that X and X 7 are 7/-related. 

On the other hand we can construct a FL/-projectable vector field 
£ on Pf which projects on X. Next we can apply Proposition 14. 161 and 
obtain the section : M f P/. Recall that X c (/3 ? (p)) = Tfe(X(p)) 
is the unique vector field on Im(/3^) which satisfies the SODE condition 
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and the equation i Xli ojl = dE L . The following lemma gives the relation 
between Im(/3g) and Qf. 

Lemma 4.17. Im(f3^) is a submanifold ofTQf. 
Proof: Since Xg verifies the SODE condition, then 
Tr Q (X 6 (p)) = r TQ (X 6 (p)) 

for any p G Im(^). 

Since X^ is tangent to Im(/3g), and since Im(^) is a submanifold of 
P f and TQ (P f ) = Q f , then Tr Q {X^p)) G TQf. 

On the other hand ttq{X^(j>)) = p G Im(^), and using the SODE 
condition we deduce that p G TQf. □ 

Remember that X$ and X are /^-related and X and X 1 are 7/- 
related, so we deduce that X^ and X 1 are /^07^-related too. Moreover, 
since Xg satisfies the SODE condition, we can find a better description 
of fa o 7^. 

Proposition 4.18. We have 

h o 7/ = 

Proof Since X^ verifies the SODE condition, then given q G Qf we 
obtain 

Tr,(X 5 (/3 c o 7/ (g))) = r TQ (X^ o 7/ (g))). 

Therefore, 

Tr Q (X ? (/3 ? o 7/ (g))) = Trg o Tfr(X(q)) = Tn(X(>y(q))) = X\q) 

where we have used that tq = n o FL and FL o = id M/ . 
On the other hand, 

r TQ (X^(3 ( o 7/ (g))) = (/% o 7/ (g)) 

Then, using the SODE condition we get X 7 = o 7^. □ 
The following corollary is immediate. 

Corollary 4.19. TTie vector fields X$ and X 1 are X 1 -related, i.e. 

X^o V {q))=TX\X\q)) 

or equivalently 

X ? (/%o 7/ (g)) = (X 7 f(X 7 (g)), 
where (X 7 ) c denotes the complete lift of the vector field X 7 . 
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4.9. Example: Lagrangian setting. 

Example 4.20. We will revisite example 14 . 1 1 1 and discuss the Hamilton- 
Jacobi problem for the Euler-Lagrange equation. The lagrangian func- 
tion is 

L{q\q\q\q 2 ) = \{q l ? + q 2 q l + q 1 q l . 
Then FL was given 

FL(q 1 ,q 2 ,q 1 ,q 2 ) = (q 1 ,q 2 ,q 1 +q 2 ,q 1 ) 
and the primary constraints are 

^ 1 {q A ,p A ) =P2~q 1 

So 

^i — {{q 1 } q 2 iP\iP'i) £ IR 4 sucn that p 2 = q 1 }- 
and we can use (q 1 , q 2 ,pi) as coordinates on Mi. 
It follows that 

El = \q l 

hi = \{pi~q 2 ) 

tu 1 = dq 1 A dp 1 + dq 2 A dq l 

Let 

h(q A ,PA) = - V 2 ) 
be an extension of the hamiltonian hi. 

It is easy to see that, at the points of Mi := lm(FL), we have 

{Q^h + u® 1 } = 

and therefore we are in presence of global dynamics 

The solution of the equation (ix^i = dhi)\ Ml is given by 

where / E C°°(Mi) 

Recall also, that a solution of our Hamilton- Jacobi problem, ^(q 1 , q 2 ) 
= (q\q 2 ,li(q\q 2 ),l2(q\q 2 )), was given by 

l{q\q 2 ) = {q\q\q\q l ) 

If we take a solution X = (pi - q 2 )^ + + + (pi - q 2 )^; we 
can compute 

d . d 
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and also 
then 

This vector field along X 7 satisfies the SODE condition. We can 
consider now the equation i^w £ = dE L 

ool = dij^dq 1 + §±dq 2 ) = ddq 1 + q^dq 1 + q'dq 2 ) 
= dq 1 A dq 1 + dq 2 A dq 1 + dq l A dq 2 = dq 1 A dq 1 . 

So, «(X7)C(X7) ojl = and dEi^X 1 ) = q 1 dq 1 {X 1 ) = and thus 

wl = dE^X 1 ). 

Therefore (X 7 ) C (X 7 ) satisfies Euler-Lagrange equations and the SODE 
condition. 
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